A first principles study of magnetism in PdsFe under pressure 
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Recent experiments on PdaFe intermetallics [Phys. Rev. Lett. 102, 237202 (2009)] have revealed 
that the system behaves like a classical invar alloy under high pressure. The experimental pressure- 
volume relation suggests an anomalous volume collapse and a substantial increase in bulk modulus 
around the pressure where invar behavior is observed. With the help of first-principles density 
functional theory based calculations, we have explored various magnetic phases (ferromagnetic, fully 
and partially disordered local moment, spin spiral) in order to understand the effect of pressure 
on magnetism. Our calculations reveal that the system does not undergo a transition from a 
ferromagnetic to a spin-disordered state as was thought to be the possible mechanism to explain the 
invar behavior of this system. We rather suggest that the anomaly in the system could possibly be 
due to the transition from a coUinear state to non-coUinear magnetic states upon the application of 
pressure. 

PACS numbers: 63.20.dk, 63.50.Gh 



I. INTRODUCTION 

The technological usefulness of the state-of-the-art 
magnetic materials depends on their properties at finite 
temperature and pressure. Fe-based alloys and inter- 
metallics have their rich history of being used as smart 
magnetic materialsi. Numerous experimental and the- 
oretical investigations have been carried out to under- 
stand the magnetic properties of pure Fe and its alloys 
and intermetallics at finite temperatures^^^ (and refer- 
ences therein). A very special characteristic of some Fe- 
based alloys is observed as the invar properties where 
the thermal expansion coefficient remains invariant with 
temperature. A wealth of literature on theoretical and 
experimental results on invar alloys exist, ^^^P Besides 
the conventional temperature invariant invar properties, 
pressure induced invar characteristics have also been 
report ed,l- 

In this paper, we will discuss an Fe-based alloy where 
pressure induced invar properties have been observed in 
experiments recently^iS, EDXAD experiments on PdaFe 
alloys have reported an anomalous volume collapse un- 
der pressure between 10 and 15 GPa^^. The EDXAD 
data were successfully fit to a equation of state which 
points towards the existence of two separate states; one 
at high pressure and the other at low pressure. Sub- 
sequent measurements at high temperatures reveal that 
at 7 GPa pressure, there is an anomalously low thermal 
expansion suggesting that the system behaves like an In- 
var alloj^ at high pressure. First principles calculations 
to understand the origin of these anomalies were rather 
inconclusive. However, in order to explain the invar be- 
havior under pressure in this system, it was conjectured 
that the application of pressure may reduce the Curie 
temperature below room temperature, thus making pos- 
sible a transition in magnetic state from ferromagnetic 
to a spin-disordered one and that the invar anomaly is a 
result of such high volume-high spin to low volume-low 



spin transition^^. Such a proposition stems from the facts 
that in Fe-rich FcPt and FePd alloys, the classical invar 
anomaly at ambient pressure was explained in terms of 
the transition from a high volume ferromagnetic state 
to a low volume spin disordered state^"^ . In case of Fe- 
based alloys with low Fe concentration, the magnetically 
ordered state was found to be stabilized and thermal in- 
var behavior was suppressed, in general. However, it was 
discovered that in case of Fe-Ni alloys with low Fe concen- 
tration, the alloy with normal thermal expansion proper- 
ties at ambient pressure exhibits invar behavior at high 
pressureJi. This explanation for this anomaly was given 
in terms of transition to non-collinear configurations at 
lower volumes due to application of pressure. On the 
other hand, Khmelevsky et alt^ investigated the mag- 
netic transitions in Feo.rPto.s under pressure and ruled 
out the possibility of making a connection between the 
observed invar anomaly and a transition to noncoUinear 
states from zero pressure coUinear one. These results, 
thus, indicate that the physical origin of anomalies re- 
lated to the magnetic properties in Fe-based alloys under 
pressure may not be unique. In this paper, we, there- 
fore, aim at understanding the magnetism in PdsFe al- 
loys under pressure by extensive quantitative analysis of 
the energetics, the magnetization, the pressure-volume 
relations with different magnetic structures. To this end 
we have performed first principles calculations with the 
coUinear magnetic structures viz. ferromagnetic, antifer- 
romagnetic, and the the spin disordered structures and 
with the spiral magnetic structures. Our results show 
that contrary to the proposition in Ref. the volume 
collapse and the subsequent invar anomaly under pres- 
sure, observed experimentally, cannot be explained by 
a transition from magnetically ordered state to a spin 
disordered state. Rather, the results indicate that the 
anomaly could be due to the transition from a coUinear 
magnetic state to noncoUinear ones. 

The paper is organized as follows. In Section II, we 
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briefly discuss the details of computational methods used. 
In Section III, we present our results on various magnetic 
structures. A detailed discussion followed by the conclu- 
sions are presented at the end. 



II. COMPUTATIONAL DETAILS 

All calculations have been performed within the stan- 
dard framework of the spin polarized version of the den- 
sity functional theor y The spin disordered magnetic 
structure has been simulated using the disordered local 
moment (DLM) formalism^^ where the magnetic disor- 
der is represented by considering the PdsFe system as a 
pseudoternary alloy Pd3Fe+Fe^_^, where x is the concen- 
tration of Fe atoms with up spin (Fe"'") and 1 — a; is with 
down spin (Fe~). The paramagnetic or the full DLM 
state, henceforth mentioned as FDLM, is obtained when 
X = 0.5. Partial disordered local moment (PDLM) is re- 
alized for X ^ 0.5. Calculations for the DLM states are 
performed by the Green's function based Exact Muffin- 
tin Orbital (EMTO)i2r2S formalism in conjunction with 
the Coherent Potential Approximation (CPA)^, a single 
site mean field formalism to treat the chemical and mag- 
netic disorder. The EMTO-CPA is a suitable method 
for the treatment of the spin disordered state as the us- 
age of conventional electronic structure methods require 
construction of prohibitively large supercell to simulate 
the magnetically disordered states. The other advantage 
of the EMTO method is that it uses two different ap- 
proaches for calculations of the one-electron states and 
the potentials:while the one electron states are calculated 
exactly for the overlapping muffin-tin potentials, the so- 
lution of the Poisson's equation can include certain shape 
approximations. An accuracy at a level comparable 
to the full-potential techniques can, thus, be sustained, 
without a significant loss of accuracy. However, since the 
present implementation of the EMTO method does not 
address noncoUinear magnetic ordering, the spiral mag- 
netic structures have been studied by the full-potential 
linearized augmented plane wave (FP-LAPW) methodSi 
using the ELK code^^ in an implementation which al- 
lows noncoUinear magnetism including spin spirals22i2i. 
In this implementation, the full magnetization density is 
used in addition to the full charge density and the full 
potential. The magnetic moment is allowed to vary both 
in magnitude and in direction inside the atomic spheres 
as well as in the interstitial regions. 

All calculations have been performed within the local 
spin density approximation (LSDA)'^^ for the exchange- 
correlation part of the Hamiltonian. For the EMTO- 
CPA calculations, the Green's functions were calculated 
for 32 complex energy points distributed exponentially 
on a semi-circular contour. The total energies were cal- 
culated by the full charge density technique^^. The s,p 
and d orbitals were included in the EMTO basis set for 
the expansion of the wave functions. The Brillouin zone 
integrations were performed with 1540 k-points. In the 
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FIG. 1. (Color online) Total energy as a function of Vred for 
FM, AFM and DLM states computed by the EMTO-CPA 
method. The energy of the ferromagnetic state is taken to be 
zero. The results of FP-LAPW calculations for FM and AFM 
states are shown for comparison. 

FP-LAPW method, the muffin tin radii of Fe and Pd were 
taken to be 1.16 and 1.27 A respectively. The plane wave 
cutoff for basis functions is RKmax = 8. The cutoff for 
the charge density was considered to be G'maa;=12. The 
number of k-points for the Brillouin zone integrations 
were 364. A broadening of 0.14 eV was used accord- 
ing to Methfessel-PaxtonSXscheme. Total energies were 
converged to less than 1 meV per atom for all the calcu- 
lations. 



III. RESULTS AND DISCUSSIONS 

In Ref. m, pressure- volume relations for PdsFe showed 
an anomalous reduction in volume from Vred ~ ~ 0-96 
to 0.91 between 10 and 15 GPa pressures, Vb being the 
equilibrium volume of the ferromagnetic state. Since 
the pressure-volume curve could be smoothly fitted to 
a Weiss like equation of state, it was proposed that the 
anomaly is related to the transition from a high volume- 
high spin state to a low volume-low spin state. On the 
other hand, the low thermal expansion under pressure 
suggested an invar anomaly and its connection to the two 
state model of Weiss-^. The authors, however, could not 
find out a stable low spin state at lower volume from first- 
principles calculations. Therefore it was proposed that 
the required low volume-low spin state could be a para- 
magnetic state. Such a transition, according to them, 
could be because of lowering of the Curie temperature 
Tc due to the application of pressure. In this section, we 
present results of extensive first-principles calculations to 
resolve these issues. 

In Fig.[Tl we present results on the relative stabilities of 
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ferromagnetic (FM), antiferromagnetic (AFM) and para- 
magnetic (PM) states by EMTO, EMTO-CPA and FP- 
LAPW methods. PM state has been modeled by a full 
disordered local moment (FDLM) treatment. The AFM 
state, considered for comparison, has Fe moments order- 
ing ferromagnetically in alternating (110) planes. We 
have done calculations for two other AFM states, in one 
of which the Fe moments order ferromagnetically in alter- 
nating (100) planes and in the other in alternating (111) 
planes. The energy of the first one, the (110) AFM was 
the lowest and hence it has been included for compari- 
son. The results of the EMTO and EMTO-CPA calcu- 
lations show that the AFM state becomes stable in com- 
parison to the FM state at Vred ~ 0.85, while the DLM 
state would be stable at a much higher compression. The 
FP-LAPW calculations predict that the AFM state be- 
comes stable at Vred = 0.9. It is to be noted that the 
first-principles projector augmented wave calculations of 
Rcf. J2i also observed stabilization of the AFM state at 
around Vred = 0.9. The proximity of the FP-LAPW and 
the projector augmented wave results are not surprising 
since both solve the Kohn-Sham equations without in- 
voking any shape approximations on the charge density 
and potential. A different Vred value for FM-AFM tran- 
sition as computed by the EMTO is due to the fact that 
certain degree of shape approximations is inherent in the 
muffin-tin orbital based methods. However, this discrep- 
ancy does not seriously affect the physics as our primary 
concern here is to find out a magnetic phase, energetically 
lower than the FM state at around the experimental pres- 
sure where the anomalous behavior in pressure-volume 
relation is found and relate the stability of the phase 
with the anomalous behavior. It is to be noted that the 
Vred values of 0.9 obtained by the FP-LAPW or of 0.85 
obtained by the EMTO for the FM-AFM transition are 
quite close to the Vred value of 0.91 where the anomalous 
volume collapse was observed in EDXD measurements. 
However, the outcome that the FDLM state is having 
a energy higher than the FM state until a much larger 
compression, rules out the proposition of Ref. [H, that 
the FM state transforms to the paramagnetic state upon 
compression with a lowering of Tc- In order to validate 
our argument, we have calculated the Tc as a function of 
the compression, as shown in Fig. [21 The has been 
estimated by the mean-field approximation as, 

2 

ksTc — - {EpDLM — Epm) (1) 

where Epdlm, Epm are the total energies of the system 
in FDLM state and in the ferromagnetic states respec- 
tively. The calculated Tc of 472K compares well with the 
experimental value of 499K^^ at zero pressure [Vred = 1)- 
The results suggest that the Tc indeed decreases with 
pressure and reduces below room temperature at around 
Vred = 0.85. This, though, may tempt one to infer that 
at this compression, the system should undergo a transi- 
tion from FM to PM state, making a correlation between 
this result and the anomalous volume collapse, observed 
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FIG. 2. Curie temperature Tc as a function of Vred computed 
by the EMTO-CPA method. 



experimentally, requires two things: first, the zero pres- 
sure volume of the PM state should be significantly lower 
than the FM state, introducing volume magnetostriction 
as is observed in the FePt and FePd alloys with high Fe 
concentrations where the invar anomaly is observed^ and 
second, the bulk modulus of the low volume state should 
be substantially higher than the high volume FM state 
as was observed experimentally^^. In Table I, we present 
results on zero pressure lattice constant and bulk modu- 
lus computed for the FM, the AFM and the DLM states. 
Results obtained by the EMTO and the FP-LAPW meth- 
ods along with the theoretical results of Ref. [13 are pre- 
sented for comparison. The results clearly show that 
the equilibrium lattice constants and the bulk modulus 
of the FM, AFM and FDLM states calculated by the 
EMTO-CPA method are nearly equal, thus leaving no 
scope of observing magnetostriction and the features in 
bulk modulus observed in experimental data of Ref. [l3- 
The FP-LAPW calculations for FM and AFM states also 
imply the same trend with regard to these two magnetic 
states. Another feature of the invar anomaly observed 
in case of FeaPt alloy was that the low volume FDLM 
state had significant reduction in the magnetic moment 
at the Fe site in comparison to that in the FM stated. 
The magnetovolume effect connected to this reduction of 
Fe moment in the FDLM state was the reason behind 
the invar anomaly^. To explore this aspect, in Fig. |31 
we show the variation of the magnetic moment at the 
Fe site with compression for the FM, the AFM and the 
FDLM states. The results suggest that the Fe moments 
are extremely localized and there is hardly any deviation 
of the Fe local moment across the magnetic orders. All 
these observations, particularly regarding the behavior 
of structural and magnetic properties in the FM and the 
FDLM states, suggest that the experimentally observed 
features under pressure in PdsFe are not due to a transi- 
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FIG. 3. (Color online) Local moments of Fe for FM, AFM 
and FDLM states as a function of compression computed by 
the EMTO-CPA method. (Inset) Total energy of LS state 
in meV per atom with respect to FM state as a function 
of compression computed by the FP-LAPW and the EMTO 
methods. 



tion from the FM to a paramagnetic (FDLM) state. Also, 
we show a comparison between the relative total energies 
obtained in EMTO and FPLAPW methods in the inset 
of Fig. [21 In both methods, FM state has energy lower 
than that of LS state and the trend is similar as a func- 
tion of compression. In Table I, we have also presented 
results on the equilibrium lattice constant and the bulk 
modulus for a low-spin(LS) state. The calculations for 
this state was done by keeping the spins collinear but by 
fixing the total magnetic moment as O.Ol/is. The consid- 
eration of such a state was part of the attempt to search 
for a low volume low spin state which may be related 
to the explanation of the anomalies observed experimen- 
tally. The results suggest that the equilibrium volume of 
this LS state is significantly lower than the FM state and 
the bulk modulus considerably higher. In Fig. |4l we plot 
the pressure-volume relations computed by both EMTO 
and FP-LAPW methods. The results show a striking 
agreement between the experiment and the theory, both 
qualitatively, as well as quantitatively when one consid- 
ers only the FM and the LS states. The pressure-volume 
relation for AFM state is also shown for comparison. It 
appears that the desired anomalous reduction in volume 
between 10-15 GPa can be explained by a transition from 
FM to LS state. However, upon calculating the total en- 
ergy of the LS state we find that the LS state is not 
stable in comparison to the FM state in the desired pres- 
sure regime as shown in Fig. |4l In fact, the energy of 
the LS state is more than 100 meV /atom higher than 
the FM state around the compression (inset of Fig. [3]) 
where the anomaly is observed and that it becomes sta- 



TABLE I. The equilibrium(zero pressure) lattice constants 
ao in A and bulk modulus Bq in GPa for different magnetic 
states computed by different electronic structure methods. 

Ref. 

do -Bq 



Magnetic 


EMTO-CPA 


FP-LAPW 


state 


ao 


Bo 


ao 


Bo 


FM 


3.888 


204.8 


3.794 


216.7 


AFM 


3.886 


202.6 


3.792 


218.8 


LS 


3.850 


229.0 


3.756 


242.3 


FDLM 


3.887 


206.1 







3.802 
3.798 
3.760 



216.8 
216.4 
242.6 



ble only at a much greater compression than the FDLM 
state. Fig. U suggests that this qualitative feature is inde- 
pendent of the method of calculation used. This is also 
reported in Ref. [H. Thus the excellent agreement be- 
tween the theory and the experimental pressure-volume 
relations can be purely fortuitous and can not be used as 
a viable explanation of the anomalies. 

In order to resolve the observed discrepancy between 
the theoretical and experimental results, we next con- 
sider the following two possibilities: (i) upon application 
of pressure, the system does not make a transition from 
FM to a FDLM state but to a partially spin disordered 
state (PDLM) and (ii) upon application of pressure, the 
system may not retain the perfect chemical order, instead 
the system may tend to be partially ordered or fully dis- 
ordered. 

To check the first possibility, we have calculated the 
total energies of the various PDLM states with rtred — 

f!: being equal to 0.6, 0.7, 0.8 and 0.9 where nK^ and 



"total 



^totai number of spin-up Fe atoms and the total 

number of Fe atoms {n^'^ -\-n^'^) respectively, being 
the number of spin-down Fe atoms. The FDLM state is 
represented by Ured = 0.5. In Fig. [SJ we present the to- 
tal energies of these spin disordered states relative to the 
total energy of the FM state as a function of compres- 
sion. The results suggest that none of the PDLM states 
is energetically lower than the FM state for any value of 
volume compression. Interestingly, one may observe that 
for Vred < 0.75, the difference in energies between the 
PDLM states with different rired values are very close in 
energy. As clear from Fig. [1] this is the volume region, 
where EMTO-CPA should yield a FM-FDLM transition. 
This suggests that the transition may not occur through 
a FDLM state but a combination of PDLM states with 
different degrees of spin-disorder. 

The second possibility as described above is explored 
by calculating total energies for partial chemical order. 
The fully ordered PdsFe has a CusAu structure with Fe 
on Au(I) sites and Pd on Cu(II) sites. In the partially 
ordered states, Pd antisites are introduced in the Fe (I) 
sites and the remaining Fe atoms are distributed on Pd 
(II) sites. In order to describe the degree of chemical 
order, the renormalized chemical long-range-order (LRO) 
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various degrees of spin disorder is not likely when one 
approaches the volume where the FM-FDLM transition 
takes place. 



10 15 20 25 
Pressure (GPa) 



FIG. 4. (Color online) Pressure-volume relations for PdsFe 
computed by the FP-LAPW method (top panel) and the 
EMTO method (bottom panel). Results of FM (blue dashed), 
AFM (black) and LS (green dashed) states are shown in the 
plot. Experimental results (red line with error bars) are in- 
cluded for comparison. 



parameter is defined in the usual way for site I as 



c{Pd)-c/{Pd) 
dPd) ' 



(2) 



where c{Pd) is total concentration of Pd in the alloy 
and c^{Pd) is a concentration of Pd atoms on Fe site 
I. Thus, S = corresponds to the fully disordered al- 
loy and S = I corresponds to the fully ordered one. In 
Figs. [5] and [71 we show the results of total energy cal- 
culations for FM and various DLM states for S = 0.5 
and S = respectively. The first noteworthy outcome of 
these calculations is that the equilibrium lattice constants 
do not change significantly with the degree of chemical 
order. The lattice constants for S — 0.5 is 3.891 A and 
that for = is 3.892 A and thus hardly differ with 
that for S — 1 given in Table I. The equilibrium bulk 
modulus for S = 0.5 is 203.9 GPa and that for = 
is 199.5 GPa, once again showing a non-significant de- 
viation from the value obtained for S — 1 (Table I). We 
observe the same trends for various DLM states as well. 
As was observed in case of 5 = 1 , no signature of volume 
magneto-striction is observed for the partially disordered 
and fully disordered states. The results shown in Figs. [B] 
and [7] suggest that same trends in relative stabilities of 
PDLM and FM states as observed in Fig.[5]are preserved. 
With increasing degree of disorder, the energy differences 
between rired — 0.8 and rired — 0.9 PDLM states for the 
compression range Vred — 0.9 — 0.85, in fact increase in 
comparison to the fully ordered state, suggesting that the 
possibility of realizing a mixed state comprising of differ- 
ent PDLM states in this compression range decreases. 
In summary, the presence of both chemical disorder and 
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FIG. 5. (Color online) Variation of the total energies for var- 
ious PDLM states as a function of compression for fully or- 
dered {S = 1) PdaFe. Different PDLM states are marked 
by the corresponding rired values. The energies are measured 
relative to the energy of the FM configuration. 




FIG. 6. (Color online) Same as in Fig. [5] but for S = 0.5, 
characterizing a partially chemical disordered state. 

These results, therefore, indicate that the degree of 
chemical order and the degree of spin disorder do not 
provide us any indication towards a transition from high- 
volume high-spin state to a low- volume low-spin state in 
PdaFc under pressure. However, the results presented so 
far were based upon coUinear magnetic configurations. It 
has been argued that the invar anomaly in Feo.64Nio.36 
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FIG. 7. (Color online) Same as in Fig. [5] but for 5 = 0, 
characterizing a completely disordered state. 



alloys is related to the transition from collinear high 
spin state to noncoUinear ones^. It is to be noted that 
the invar anomaly accompanied by the anomaly in the 
pressure-volume relation and in the variation of bulk 
modulus has been experimentally observed for both am- 
bient and high pressures for various compositions, includ- 
ing compositions with low Fe concentrationsii. 




[qOO] 



FIG. 8. Total energy as a function of spiral vector q in units 
of 27r/a along [qOO] direction for different compressions. In 
the parentheses, the values of K-ed are shown. 

Following this idea, we study the spiral magnetic states 
in PdsFe. Specifically, we look for the possibility of non- 
coUinear ordering by studying the energetics of spiral con- 
figurations. The total energy of the system is calculated 
as a function of the spiral wave vector q for different 
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FIG. 9. Same as in Fig. [8] but for q along [qqO] direction. 
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FIG. 10. Same as in Fig. [8]but for q along [qqq] direction. 



values of the compression K-ed- In Figs. 8, 9 and 10, 
we present the results for the energetics along the high- 
symmetry directions [100], [110] and [111] respectively 
for different compressions. Polesya et al.'^^ have reported 
the spin spiral calculations for both ordered and disor- 
dered FdsFe under ambient conditions. Our calculated 
results along [100] direction for Vred = 1-00 are very sim- 
ilar to their data. The behavior of the spin spiral curve 
is explained in terms of the induced moment on Pd sites. 
Also, the reduction of Pd moment along the [100] direc- 
tion (data not shown here) is similar in both studies. 

The analysis of the [100] and [110] spin spirals show 
that along the [110] direction, there is a tendency to- 
wards the stabilization of the [0.5 0.5 0] spin spiral with 
respect to the [000] one starting from Vred = 0.96, mak- 
ing it lower in energy than [000] one at Vred = 0.92. This 
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is consistent with the results of the colhnear calculations 
shown in Fig. [T] where the AFM state becomes energeti- 
cally stable than the FM around Vred = 0.9. Although, 
no incommensurate spin wave was found to be the low- 
est in energy, the energy of the [0.4 0] spin spiral also 
started to decrease significantly from Vred — 0.96, finally 
producing a total energy lower than the FM phase(q=0 
0) for Vred — 0.9. Further compression resulted in the 
complete stabilization of various incommensurate phases 
over the FM phase. Analysis for [q 0] spin spirals sug- 
gest that although no incommensurate phase or the AFM 
phase ((7=0.5) is lower in energy than the FM phase (g=0) 
upto Vred — 0.92, the energy of the spirals for q > 0.4 
starts decreasing continuously and at Vred = 0.9, they are 
within 4 meV per atom in energy to the FM state. For 
[111] direction, we do not see any sign of stabilization of 
any incommensurate or AFM phase (q=0.5 0.5 0.5) with 
respect to FM phase. However, the energies for q < 0.2, 
are less than 5 meV per atom with respect to that in the 
FM phase. 

The following picture emerges from the assimilation 
of all the results on collinear magnetic states and spin 
spirals: first, the anomaly observed in PdaFe under pres- 
sure is not a simple effect of transition from high volume 
FM state to a low volume paramagnetic or any low spin 
collinear state, and second, for the range of compression, 
where the anomaly is seen, there is a close competition 
between the FM, the AFM and other incommensurate 
states. This can be interpreted from the results presented 
in Fig.[T]and Figs. 8-10. This close competition points to 
the fact that the tendency of the system is to drift from 
a collinear high spin state to increasingly noncollinear 
states as a function of increasing pressure. The anomaly 
in the pressure-volume relation and the anomaly in the 
bulk modulus for different compressions, observed exper- 
imentally, act as a confirmation to this interpretation. 
This similar anomaly observed in Feo.65Nio.35 at ambient 
prcssurc^and for Feo.20Nio.8O sA, finite pressures^^ was ex- 
plained by the transition from FM to noncollinear states. 



IV. CONCLUSIONS 

We have performed extensive first-principles calcula- 
tions on different magnetic phases of PdsFe under ex- 
ternal pressure to understand the invar related pressure- 
volume anomaly in a recent experiment on PdsFe. The 
results are analyzed on the basis of relative stability of 
various magnetic phases as a function of compression. 
Contrary to the conjecture made in Ref. [H, we find 
that the experimentally observed invar anomaly in the 
pressure- volume relations and bulk modulus, can not be 
explained based upon either the destabilization of the 
equilibrium FM phase to a collinear low-spin phase or to 
a paramagnetic or partially spin disordered state under 
application of pressure. Non-collinear calculations with 
spin spirals also could not resolve the issue by stabiliza- 
tion of any incommensurate excitation under compres- 
sion. However, our results suggest that in the pressure 
range where the anomaly is observed in experiments, the 
FM, AFM and some incommensurate magnetic phases 
compete closely, indicating that the collinear magnetic 
structure may be destabilized under external pressure, 
making the transition to non-collinear structures by the 
anomalous reduction of volume and increase in bulk mod- 
ulus as has been observed earlier in other systems. A final 
resolution on this issue can be made by first-principles 
calculations on complete non-collinear magnetic struc- 
tures as a function of pressure. This will be addressed in 
a future communication. 
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